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An overdoped YBa2Cu307 single crystal was studied by SQUID and torque magnetometry in 
order to investigate the temperature dependence of the anisotropy parameter close to the transition 
temperature T c (0.87 T c < T < T c ). Angle dependent torque measurements were performed and 
analyzed with the widely used Kogan model [Phys. Rev. B 38, 7049 (1988)] as well as with an 
extended model by Hao and Clem [Phys. Rev. Lett. 67, 2371 (1991)], taking into account the 
influence of the vortex cores on the magnetization. Both approaches yield similar results, with an 
out-of-plane anisotropy parameter around 6.5 which slightly increases with decreasing temperature, 
and a temperature independent in-plane anisotropy parameter "/ a t = 1.12(5). 
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I. INTRODUCTION 

Since the discovery of high-temperature superconduc- 
tivity in the cuprates P the anisotropic magnetic prop- 
erties of layered superconductors were extensively stud- 
ied (see e.g. Refs. [^HS]). All superconductors with a 
transition temperature T c > 30 K have a layered struc- 
ture. In particular, the superconducting gap of cuprates 
was found to be strongly anisotropic due to the crystal 
structure consisting of weakly coupled superconducting 
Cu02 planesP^ It is interesting to investigate how the 
anisotropic properties change as a function of thermody- 
namic parameters and doping within a particular family 
of cuprates, and to explore their common features by 
comparing various families. 

The gap structure can be probed directly by excit- 
ing superconducting carriers. Importantly, the energy 
needed for this, i.e. the energy gap, may be accessed by 
probing the magnetic penetration depth AP^In a layered 
superconductor the gap structure is strongly anisotropic, 
thus A is anisotropic as well. The magnetic penetration 
depth related to a supercurrent flowing along the i-axis 
(i = a,b,c) is denoted as Aj, and the penetration depth 
anisotropy between two crystallographic directions i and 
j is jij = Xi/Xj. 

In the anisotropic Ginzburg-Landau theory, which is 
the most commonly applied phenomenological descrip- 
tion of layered superconductors, the anisotropy is de- 
scribed by the temperature independent effective mass 
anisotropy (assuming XJXj = (m^/m*) 1 / 2 = H^/H^, 

where ifjjf is the upper critical field along the fc-axis). 
However, a temperature dependent anisotropy was ob- 
served in various superconductors, especially in the two- 
gap superconductor MgB 2 (see Refs. [7] and Q2]), and 
was explained as a consequence of the presence of two 
superconducting gaps. A similar temperature depen- 
dence was also observed in iron-based superconductors,^ 
in which evidence for two-band superconductivity was 



provided by seve ral e xperiments, i nclud ing point con- 
tact spectroscopy^^ and ARPESpEU Multi- gap su- 
perconductivity seems to be more common than first 

expected, as indications of it were also observed in 

cupratesPSEni ft m ay be related to the temperature de- 
pendence of the anisotropy} 2 ^ as in the case of MgB^and 
the iron-based superconductors.^ However, there may 
be other reasons for t his t emperature dependence: the 
anisotropy of the gapj 11 * 22 ! the anisotropy of the Fermi 
surface, 23 or strong coupling.- ' A temperature depen- 
dent anisotropy parameter was also observed in cuprates 
(see e.g. Refs. M IT51 [251 and 125)) . This rises the question 
whether the temperature dependence of the anisotropy is 
a common property of all layered high-T c superconduc- 
tors, and how it is linked to the gap structure. 

A recent study of the cupratc superconductor 
SmBa2Cu307_ ( 5 facing this question was limited to the 
underdoped region onlyP It was shown that the temper- 
ature dependence of the anisotropy is more pronounced 
for samples with lower oxygen content. Such samples 
are characterized by a well developed pseudogap, i.e. an 
additional energy scale which may play a similar role in 
the development of the temperature dependence of the 
anisotropy as the multi-gap structure in MgB2 and iron- 
based superconductors. Therefore, it is very important to 
perform reliable studies of the temperature dependence 
of the penetration depth anisotropy for optimally doped 
and overdoped cuprates. In this doping range the pseu- 
dogap vanishes or eventually overlaps with the supercon- 
ducting gap. 

Taking all of the above into account, we decided to 
study the temperature dependence of the anisotropy of 
a detwinned, almost fully oxygenated, overdoped sin- 
gle crystal of YBa2Cu30*7_,5. This system exhibits an 
anisotropic energy gapp3 and several experime nts in di- 
cate an order parameter of s+d wave symmetry! 18 ! 28 ! 

Here, we report on torque measurements of the 
anisotropy parameter of an overdoped YBa2Cu307 single 
crystal. Torque magnetometry provides a direct method 
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to study the anisotropic magnetic properties of supercon- 
ductors, contrary to methods measuring physical quan- 
tities separately along different crystallographic direc- 
tions from which the anisotropy parameter is determined. 
An analytical approach for the analysis of experimen- 
tal data based on the solution of the Hao-ClenP2l func- 
tional is applied, which allows to investigate anisotropic 
extreme type-II superconductors beyond the London ap- 
proximation. For simplicity, the London approxim ation 
of the anisotropic Ginzburg-Landau theory (AGLT)j 30 * 31 ! 
in which simplifications of the geometry of the vortex 
structure are made, is often used for analyzing experi- 
mental data. However, as discussed by Hao and CleirPH 
this approximation may not necessarily be adapted to the 
interpretation of magnetization measurements, and thus 
both approaches are compared in this work. 

Section II gives a brief review of the London and of the 
Hao-Clem models in connection with the torque magne- 
tometry technique used in this work. The experimental 
details are described in Sec. III. The results and the dis- 
cussion are presented in Sec. IV, followed by the conclu- 
sions in Sec. V. 



II. LONDON AND HAO-CLEM MODELS 

The angular dependent magnetization ivl of a sample 
with volume V and magnetic moment is derived from 
the free energy F of an anisotropic superconductor in the 
mixed stated 

-1*.* (!) 



The magnetic torque 



T = -V ff F = 



M x H 



(2) 



is related to M and the angle 9 between and the crys- 
tallographic c-axis. 

A direct calculation of F within AGLT is not trivial, 
since F depends on the exact distribution of vortices and 
thus on the local magnetic induction B(H) inside the 
superconductor. However, F can be expressed within 
the so-called London limit, assuming that the influence 
of the finite vortex core size can be neglected. This is 
valid if the vortex core size is very small compared to 
the vortex itself, i.e. the penetration depth is much 
larger than the coherence length. Anisotropic supercon- 
ductors exhibit disctinct magnetic properties along the 
principal axes a, b, and c. In layered superconductors the 
largest anisotropy is observed between the c-axis and the 
layers (afe-plane). Therefore we may approximate the 
orthorhombic structure of YBa 2 Cu 3 07 by a tetragonal 
one j33| introducing the anisotropy parameter 
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where the in-plane magnetic penetration depth A a h = 
v 7 A a Af>. The magnetization M and the torque r are de- 
rived in the so-called Kogan modeP^l as 



Mi, (6, H) = 



$oe(0) 



hi 



8n^X 2 ab \e(9)H 



(4) 
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,H) = - 



V<S> a H 
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Y 



Bin(2fl) 



In 



e(9)H 



(5) 

Here, the index L indicates the London approach, <&o 
is the magnetic flux quantum, and e(9) is the angular 
scaling function 



f 



e{9) = J cos 2 (9) + ^ sin 2 (9). 



V 



(6) 



The parameter rj accounts for uncertainties due to the ap- 
proximation of the London limit (e.g. the neglected sup- 
pression of the order parameter inside the vortex cores). 

Hao and ClerrP^ showed by analyzing the free energy 
within AGLT that the parameter rj cannot be constant in 
the entire magnetic field range H c % < H < H C 2, which is 
also evident from more recent theoretical worlfSl beyond 
the Hao-Clem mo del. Th e correct functional form of F 
by Hao and ClenpM5136| incorporates in the expression 
for M and r the empirical functions a(h) and /3(h), where 
h denotes the reduced field 



h(9) = 



H 



H c2 {9Y 



(7) 



Their generalized treatment of the mixed state of a 
superconductor, which includes the vortex core contribu- 
tion to the free energy functional, yields a more realistic 
formula for the magnetization^ 



M w (0,H) = -a(h(9)) 



8^ Xl b { h{9) ) 



(8) 



where the index HC indicates the Hao-Clem model. Ac- 
cording to Eq. ([2]), the torque is written as 



mc(0,H) = -a(h(d)) 



V<S> H 



(9) 



e(9) \ h{9) ) 



Here, taking into account the suppression of the or- 
der parameter in the vortex core leads to a modification 
of Eqs. (ij) and ^ by including the functions a(h) and 
13(h). These functions account for the correction of the 
in-plane magnetic penetration depth \ a b and the c-axis 
upper critical field H^, respectively. For a(h) — 1 and 
(3(h) = 77, Eqs. (|]) and @ reduce to Eqs. Q and ^ 
of the London limit. Within the HC treatment no ana- 
lytical formulas for a(h) and P(h) can be derived easily. 



3 



However, for a Ginzburg-Landau parameter k » 1 the 
following values for a and (3 are founcP*! 



0.02 < h < 0.1 : ~ 0.84, /3(h) ~ 1.08 (10) 
0.1 < ft < 0.3 : a(ft) ~ 0.70, /3f» ~ 1.74 . (11) 

It is clear that although a and /3 are assumed to be 
constant in the London limit, they are field dependent 
and may vary considerably with magnetic field (a and (3 
are fully determined by the reduced field ft)p2l 

In Fig. [I] we present the numerically calculated field 
dependence of the reduced magnetization Mnc(h) / H C 2 
and compare it with the empirical Eq. Q in order to 
extract a(h) and /3(h). The quantity kMhc (h)/H C 2 
for 2 < k < 200 is presented in Fig. flk (for clarity, 
kMhc (h)/H C 2 is shown instead of M^q(Ji) / H c2 ). Obvi- 
ously, Mnc(ft) / Hc2 strongly depends on k. The functions 
a(h) and /3(h) are presented in panels b) and c), respec- 
tively. For k > 50, the functions a(h) and (3(h) become 
essentially independent of k. The derived a(h) and j3(h) 
are in good agreement wi th t he values estimated by Hao 
and CleirP given in Eq. ilOl. 



Analyzing magnetic torque experiments by means of 
the above described theoretical model by Hao and Clem, 
one should note that the parameter k is the isotropic 
Ginzburg-Landau parameter k = A/£, where £ is the co- 
herence length, and A is the magnetic penetration depth. 
However, for a layered superconductor, the Ginzburg- 
Landau parameter is anisotropic: k has to be replaced by 
n(9) = K c /e(9) 2 , where k c = A c /£ c = j 2 \ a b/£ab = J 2 Kab- 
The functions a(h) and (3(h) then depend on the angle 9 
not only via h(9), but also via k(9). For YBa2Cu307 in 
a field of 1.4 T at T — 80 K, h(9) varies approximately 
between 0.1 (0 = 0°) and 0.01 (9 = 180°) when the field 
is turned from the c-axis to the a6-plane. In this case, a 
and (3 strongly depend on h(9) (see Figs. [T]b and c), but 
for k(9) > 50 they are independent of n(9). Using the 
values for X ab (T = 0) and £ ab (T = 0) for YBa 2 Cu 3 7 
from Refs. I3T1 and I551 one gets K ab ~ 400, which means 
that k(9) > 400 for all 9. Thus, for YBa 2 Cu 3 07 the pa- 
rameter k(9) has negligible influence on a and (3. As a 
consequence, taking into account the anisotropy on the 
Ginzburg-Landau parameter k in the HC model does not 
lead to a more reliable determination of a and (3. 



III. EXPERIMENTAL DETAILS 

The sample studied is an overdoped detwinned sin- 
gle crystal YBa2Cu307 grown in BaZr03 crucibles, with 
dimensions 130x160x50 /Ltm 3 and T c ~ 88 K. Crystal 
growth in BaZr03 yields samples of highest puritypS To 
fully oxygenate the crystal a high pressure annealing was 
performed at 300 °C in 100 bar of oxygen. Samples 
produced in this way show no anomalies, e.g. fishtail 
effect pQ^nd have very low pinning A Quantum Design 
MPMS XL SQUID magnetometer was used to determine 
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FIG. 1. (Color online) a) kMhc(/i)/-I?c2 for various values 
of k calculated using the HC model, b) Empirical parameter 
a(h) extracted from Mnc(h)/ H C 2- c) Empirical parameter 
(3(h) extracted in the same way as a(h) in panel b). Both 
a(h) and /3(h) are essentially independent of k for k > 50. 



T c (Fig. [2]). The temperature dependence of the magnetic 
moment m was measured in a small field fioH = 1 mT 
parallel to the a6-plane in zero field cooled (ZFC) and 
field cooled (FC) mode. The small difference between 
m(T) obtained in the two modes and the sharp transi- 
tion observed indicate a good quality of the crystal. 

The torque measurements were carried out using a 
home-made torque magnetometerP The piezoresistive 
sensor used consists of a platform connected to piezore- 
sistive legs which are bent when the sample mounted on 
the platform undergoes a torque. The resulting resis- 
tance change in the piezoresistors is detected by a Wheat- 
stone bridge. The read-out voltage is proportional to the 
torque magnitude r. The small dimensions needed for 
the sample allow the study of high-quality single crys- 
tals. 



IV. RESULTS AND DISCUSSION 

The torque measurements were performed in an exter- 
nal field of 1.4 T and 1 T in order to check for a possi- 
ble field dependence of the anisotropy. Since the sample 
has diamagnetic and anisotropic properties, its magne- 
tization M is not quite aligned with the field if, which 
results in a torque ~? oc M x {jl 1% according to Eq. (Jij). 

In general, the torque signal is distorted by pinning ef- 
fects: the vortex cores are pinned by defects in the sam- 
ple, in which superconductivity is more easily suppressed. 
Consequently, the sample is not at thermodynamic equi- 
librium during the time span of one measurement. As a 
result, the torque signals are different for angular fields 
measurements in opposite directions. In order to get re- 
versible angular dependent torque data, a "vortex shak- 
ing" technique^ was used. In this technique a small AC 
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FIG. 2. (Color online) Magnetic moment m(T) of the 
YBa2Cu307 single crystal, measured in a magnetic field 
fioH = 1 mT parallel to the afe-plane. In the ZFC mode 
the field is applied once the sample is cold, whereas in the FC 
mode the sample is cooled while the field is applied. The ZFC 
and FC magnetization curves show a sharp transition with a 
transition temperature T c = 88 K, indicating a high quality 
of the crystal. 
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FIG. 3. (Color online) Magnetic torque r as a function of 9 for 
single crystal YBa2Cus07 with and without vortex shaking 
at 83 K and 1.4 T. The shaking removes the irreversibility be- 
tween the increasing angle (up) and decreasing angle (down) 
measurements. The Kogan and the HC models both describe 
the data equally well. 



field ( ~ 200 Hz, 1 mT) is applied perpendicular to the 
main external field in order to shake the vortices out of 
their pinning sites (Fig. |3| . 

The temperature range of the angular measurements 
was 77 K to 86 K. The lower temperature bound was 
chosen such as to avoid the lock-in effect (also known as 
intrinsic pinning) J^2l which influences the torque in a way 
not accounted for in Eq. pi. When the external field 
direction becomes close to the etfr-plane (9 = 90°), the 
magnetization abruptly aligns with the a6-planes in order 
to minimize the magnetic energy in the superconducting 
state. However, it "jumps" back outside the planes when 
the external field direction is sufficiently away from the 
afo-plane. The upper temperature bound was chosen such 
as to avoid fluctuation effects^ close to T c . Fluctuation 
effects are not taken into account in the mean-field ap- 
proximation of the models considered here. 

YBa2Cu307 has an orthorhombic structure. Taking 
into account that in this case A ^ A&, one has to replace 
7 by 7ca or 7cb in Eqs. @, and @, with the 

magnetic field direction in the ac- or 6c-planej 4 l 43 l and 
A a 6 = VAqAj, is not equal to A a and A& as in the tetrago- 
nal case. In order to check the validity of this tetragonal 
approximation, measurements as a function of angle were 
performed in both the &c-plane and ac- plane (see Fig. [4|. 
As expected, the data are similar for both orientations, 
thus allowing the analysis within a tetragonal model. 

The torque data were analyzed with the HC and 
with the Kogan model. In order to reduce the num- 
ber of free fit parameters, the upper critical field was 
fixed in the fitting procedure according to a Wcrthamcr- 
Helfand-Hohcnberg (WHH) temperature dependence^ 

suitable for YBa 2 Cu 3 7: 38 HoH& - -1-9 T/K • (T — T c ). 
The paramagnetic background signal x(FiJ 2 /2)sin(2#) 
present in the torque data was subtracted using the 
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FIG. 4. (Color online) Angular dependent torque measure- 
ments of a single crystal YBa2Cu307 taken in the tempera- 
ture range between 77 K and T c at 1.4 T. (For clarity not all 
temperatures are shown; the arrows indicate increasing tem- 
perature). These raw data include a sinusoidal background, 
a) Measurements with H in the 6c-plane. b) Measurements 
with H in the ac-plane. 



method described in Refs. [45] and |46l Figure [5] shows 
the temperature dependence of the anisotropy parame- 
ter 7 c b as determined from the torque data using the two 
models. As evident in Fig. [5^i, both models yield very 
similar values for "f c b (within 2% accuracy). Moreover, 
the results depend only weakly on the value taken for 
H^2 (Fig. [5)3) and on the external field (Figs. [5^ andjijl). 
The errors of the fit parameters 7^ and X a b were esti- 
mated with a Monte-Carlo method: different fits were 
performed for randomly sampled points within the ex- 
perimental error of the measured data points. The final 
values of the parameters 7 and A were taken as the av- 
erage values obtained by this procedure, and their errors 
were defined as twice the standard deviation of these re- 
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FIG. 5. (Color online) Overview of the results of the angu- 
lar dependent magnetic torque measurements of single crys- 
tal YBa2Cu307 for various temperatures and fields, using the 
models described in the text, a) Comparison of 7 C (, for the 
HC and Kogan model at fixed hqH = 1 T and fixed upper 
critical field (fiodH^ / dT = -1.9 T/K). The Kogan model 
yields a slightly larger anisotropy parameter than the HC 
model, b) -y c b for the HC model at hqH = 1 T and vari- 
ous upper critical fields. The parameter ftodH^ /dT does not 
change the shape of 7(T). c) jd, for the Kogan model at 
(j,qH = 1 T and noH = 1.4 T and fixed upper critical field 
(HodH^/dT = -1.9 T/K). d) -y cb for the HC model, same 
conditions as in panel c). Panels c) and d) show that the field 
dependence of 7 is only marginal. 



suits. The estimated error bars are smaller than the size 
of the data points. 

Since the anisotropy parameter is only weakly field de- 
pendent (see Figs. [5J; and [5ji), we take as the final 7 
value the average for 1.4 T and 1 T.The corresponding 
temperature dependences of 7 CQ and 7^ are shown in 
Fig. [6^,. The free fit parameters in Eqs. ([5| and ^ are 
the anisotropy parameter 7^ and the in-plane magnetic 
penetration depth A a ;,. Since the volume V of the sample 
is not known precisely, the extracted value for A a f, may 
deviate from the intrinsic value. However, the shape of 
A a b(T) reflects the true temperature dependence, because 
the volume V is only slightly temperature dependent. 
Figure |6j; shows the temperature dependence of 1/A^ b 
as estimated from the torque data using the Kogan and 
HC model. Over the temperature range studied, 7 ca as 
well as 7 c h slightly increase with decreasing temperature, 
whereas the in-plane anisotropy parameter 7^ is tem- 
perature independent, in fair agreement with previous 
/iSR measurements of the magnetic penetration depth 
obtained for a similar sample.- 8 However, since ^ a b — 1 
one should note that it is difficult to draw definite conclu- 
sions about its temperature dependence. The tempera- 



ture dependence of the magnetic penetration depth from 
the previous ^SR studjffHl was measured along the three 
principal crystallographic axes and was interpreted in 
terms of a mixed order parameter of s+d wave symmetry. 
The values of j ca , ~/ cb - and j a b determined at T ~ 80 K, 
together with the values obtained by various experimen- 
tal techniques at different temperatures are summarized 
in Table [T] The small differences in the values may be 
due to the different temperature ranges, the experimental 
techniques used, or slight differences in the doping of the 
samples. A determination of the anisotropy from the ra- 
tio of the penetration depths requires a careful evaluation 
of A c and A a f,. Any misalignement of the sample with the 
applied magnetic field will result in an underestimation 
of A c and the deduced 7. In torque measurements, how- 
ever, the anisotropy is extracted from a fit to the data, 
without orientation issues since the model describes the 
variation of torque in the full angular range. Therefore, 
the obtained values are much more reliable. The vortex 
shaking technique allows us to avoid an overestimation of 
the anisotropy due to pinningP3 The parameters 7 ca and 
7 C 6 are slightly different, because of the orthorhombicity 
of the crystal structure. The torque data were analyzed 
here under the assumption that the field and penetration 
depth anisotropy parameters are equal. It is possible to 
generalize this analysis to the multi-gap case, where these 
parameters are not equal.^ However, such an analysis of 
the present torque data would not provide reliable results 
due to the too large number of fit parameters. 



TABLE I. Comparison of anisotropy parameters of 
YBa2Cu307_i obtained by various experimental techniques 
at different temperatures. 



technique 



T (K) 7aft 



fob 



7c 



Ref. 



low-energy (iSR 

SAN^l 

AtSR 

MSR 

specific heat 
torque (shaken) 
torque 
torque 



1.19(1) 

1.5 1.18(2) 

10 1.15(21 4 - 2 ( 5 !^ 

80 l.l(lf] 4.5(lf 

70-9(f] . 5.3(5, 

80 1.12(5) 7.00(5) 

90 1.18(14) 8.95(76) 7.55(63) 

93 1.1(2) 7.3(5) 6.6(5) 




[37] 

48 

M 
M 

\m 

this work 

43 



a small angle neutron scattering 
b estimated from Ref. 1181 

c temperature not specified; out-of-plane anisotropy parameter 
determined from H c 2 



Although no temperature dependence of the out-of- 
plane anisotropy parameter for Pr-doped YBa2Cu30 , 7_ ( 5 
was foundP it was noted that such a dependence 
cannot be ruled out due to the narrow temperature 
range studied (82 - 88 K). A temperature indepen- 
dent out-of-plane anisotropy parameter was also ob- 
served for HgBa 2 Ca3Cu40iop' where the studied tem- 
perature range was very narrow as well. In contrast, a 
pronounced temperature dependence of the out-of-plane 
anisotropy parameter was seen in MgB 2 , 7 which was 
consistently described in the framework of two-gap su- 
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FIG. 6. (Color online) a) Temperature dependence of 
the anisotropy parameters y ca and j c b for single crystal 
YBa2Cu307, obtained by averaging 7 for both measured fields 
(1 T and 1.4 T). b) Temperature dependence of the in- plane 
anisotropy parameter y a b. The dotted line is the average 
j a b = 1.12(5). c) Temperature dependence of 1 / A^ b for mea- 
surements with H parallel to the ac- and be- planes. 

perconductivity. The temperature dependence of the 
out-of-plane anisotropy parameter found later in iron- 
based supercondu ctors was also attributed to multi-gap 
superconductivity!^^ This may suggest that the present 
results are a signature of two-gap superconductivity in 
YBa 2 Cu307, as previously proposed in Ref. IS How- 
ever, we note that the temperature dependence of the 
out-of-plane anisotropy parameter observed for cuprates 



is extremely sensitive to the oxygen content. A well 
pronounced temperature dependence of the anisotropy 
for strongly underdoped sampled becomes very weak for 
overdoped YBa2Cu307. This may be related to the evo- 
lution of the pseudogap with doping in YBa2Cu307_a. 
Moreover, it suggests that an additional energy scale 
to the superconducting energy gap in the system is 
necessary to get a temperature dependent out-of-plane 
anisotropy parameter in layered superconductors. Such 
an energy scale may originate from the multi-gap nature 
of superconductivity in MgB 2 and in pnictides and from 
the appearance of the pseudogap in cuprates. 

V. CONCLUSIONS 

The magnetic torque of an overdoped YBa2Cu3C>7 sin- 
gle crystal was investigated at temperatures close to T c 
in magnetic fields of 1 T and 1.4 T. In the temperature 
range 0.87 T c < T < T c , the anisotropy parameters j ca 
and jcb were found to increase by more than 10% with 
decreasing temperature, but no field dependence was ob- 
served. In contrast, the in-plane anisotropy parameter 
"f a b exhibits no temperature nor field dependence. The 
values of j ca , j c b and 7^ are in good agreement with 
those reported previously (see Table [I]) . The analysis of 
the torque data with the Hao-Clem model yields within 
2% the same results as the simpler Kogan model. The 
Hao-Clem model does not provide new information on 
the vortex state of YBa2Cu3C>7 in the present study. 

The weak temperature dependence of the out-of-plane 
anisotropy parameter may indicate the presence of two 
energy scales in the superconducting behavior, related to 
multi-gap superconductivity or one-gap superconductiv- 
ity with a pseudogap. To clarify this hypothesis more 
experimental work is required. 
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